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Abstract 

We discuss the large scale effective potential for elastic objects (manifolds) in the presence of a 
random pinning potential, from the point of view of the Functional Renormalisation Group (FRG) 
and of the replica method. Both approaches suggest that the energy landscape at large scales is a 
succession of parabolic wells of random depth, matching on singular points where the effective force 
is discontinuous. These parabolas are themselves subdivided into smaller parabolas, corresponding 
to the motion of smaller length scales, in a hierarchical manner. Consequences for the dynamics of 
these pinned objects are underlined. 
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1 Introduction 



The physics of elastic objects pinned by random impurities is certainly one of the most topical current 
themes of statistical mechanics. The problem is of fundamental importance both from a theoretical 
point of view (many of the specific difficulties common to disordered systems are at stake) and for 
applications: the pinning of flux lines in superconductors jl], ^, of dislocations, of domain walls in 
magnets, or of charge density waves |5|, controls in a crucial way the properties of these materials. 
Interestingly, this problem is also intimately connected to surface and crack growth ||^] and to 
turbulence 

Two different general approaches have been proposed to describe the statics of these pinned mani- 
folds, for which perturbation theory badly fails. The first one is the 'functional renormalisation group' 
(FRG) which aims at constructing the correlation function for the effective pinning potential acting on 
long wavelengths using renormalisation group (RG) ideas ^, ^ . The second is the variational replica 
method which combines a Gaussian trial Hamiltonian with 'replica symmetry breaking' to obtain re- 
sults in the low temperature, strongly pinned phase [0, 12, ITs). Although many of the results of these 



two approaches actually turn out to be similar [14, 15, 11, 16, 17, 13|, the feeling that the link between 
them is missing is rather widespread, reflecting the fact that our present general understanding of 
disordered system is still incomplete. 

The aim of this letter is to unveil precise connections between these two (sometimes presented as 
conflicting |l^, |l^) theories. We show that both formalisms are indeed struggling to describe an 
awkward reality: the effective, long wavelength pinning potential has the shape drawn in Fig.l. It is 
a succession of parabolic wells of random depth, matching on singular points where the effective force 
(i.e. the derivative of the potential) is discontinuous. These discontinuities induce a singularity in 
the effective potential correlation function, and are encoded in the replica language by the RSB. The 
replica calculation furthermore provides an explicit construction of this effective (random) potential, 
and hence, in turn, information on the statistics of - say - the depth of the potential minima. The 
replica calculation might also shed light on the domain of validity of the FRG, by making more explicit 
the assumptions on which the latter relies. 

Apart from the satisfying possibility of reconciling two rather different microscopic methods, we 
believe that our construction is very useful to understand the dynamics of such objects. For example, 
their relaxation can be analyzed in terms of hops between the different minima ('traps'), corresponding 
to metastable long wavelength configurations. The statistics of barrier heights control the trapping 
time distribution, and hence the low frequency response and its possible aging behaviour []l9| , pO| ]. 
Another interesting situation is the zero temperature depinning transition induced by an external 
driving field, which has recently been investigated, again using RG ideas for expanding around a 



mean-field limit |21, 22, 23|. However, the results depend on the form of the pinning potential in this 
mean-field limit. The correct form was surmised by Narayan and Fisher |2^ to be the 'scalloped' 
potential of Fig 1. Our calculation, to some extent, confirms their intuition. 

The model we consider is the (by now standard) Hamiltonian describing pinned elastic manifolds: 



+ Fo(x,</>(x)) 



(1) 



where x is a D-dimensional vector labelling the internal coordinates of the object, and i;^(x) an A^- 
dimensional vector giving the position in physical space of the point labelled x. Various values of D 
and N actually correspond to interesting physical situations. For example, D = 3, N = 2 describes the 
elastic deformation of a vortex lattice (after a suitable anisotropic generalisation of Eq. (1)), D = 2, 
N = 1 describes the problem of domain walls pinned by impurities in 3 dimensional space, while D = 1 
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corresponds to the weU-known directed polymer (or single flux line) in a + 1 dimensional space. 
The elastic modulus c measures the difficulty of distorting the structure, and Vb(x, (^(x)) is a random 
pinning potential, which we shall choose to be Gaussian with a short range correlation function: 

Voi^Jm^'J'), = NW6^{^ - x')i?o ' (2) 

where W measures the strength of the pinning potential. In the following, we shall choose for conve- 
nience Roiu) = e^p(~2A^) where A is the correlation length of the random potential. 

One aim of the theory is to understand how the microscopic pinning potential will affect the elastic 
manifold on long length scales, relevant for macroscopic measurements. In other words, one would 
like to construct the effective pinning potential seen by a low wavevector mode of the structure, after 
thermalizing the modes with shorter length scales. Both the FRG and the replica approach propose 
an approximate construction of this effective potential which we now discuss and relate. 



2 The Functional Renormalisation Group. 

In the spirit of the momentum shell renormalisation group, the FRG method consists in writing down 
a recursion relation for the correlation function of the potential acting on 'slow' modes 0<, after 'fast' 
modes ^> (corresponding to wavevectors in the high- momentum shell [A/6, A]) have been integrated 
out using perturbation theory. This procedure has been addressed in considerable detail in Ref. |1C], we 
present only a brief description of the calculations. At zero temperature the renormalized Hamiltonian 
is defined by Hr[^^] = ^ ^|V(?<|2 + Vk[$<] and 

Vr[<A<] =min£|i|V0>|2 + y(<^< + 0>,x)|, (3) 

where the original field (p = </)<+</)> has been split into low (</><) and high ((/>>) momentum components. 
The renormalized Hamiltonian H^^ thus describes the long-distance physics of modes with momenta 
k < A/6, where A is the original short-scale cutoff, and the rescaling factor 6 > 1. The FRG proceeds 
to determine the minimum in Eq. (|3|) perturbatively in cj)^. The extremal condition may be expanded 
in 4>y as 

-vV> = -a,y(0< + <?>,x) 

« -a,y(0<,x)-a,a,y(<?<,x),/.i. (4) 

where di = Defining the Fourier transform = didj ■ ■ ■ J^V{(j)<^,x)e~^^'^, the approximate 

solution is 

AV>,k - -V^: + / v^^v^. (5) 



Inserting this solution into the energy (Eq. (|3D) gives 

where is restricted to the high-momentum shell. If 4><c{^) is constant over regions of size £, 
this can be rewritten as an integral of a local potential, up to small errors of order V/j ((/><) — 
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/ dx V/j ((/><, x). Thus, in the long wavelength limit, the renormalized Hamiltonian is well-described 
simply by a renormalized potential. Its connected correlations can be calculated from the expression 



/yi2 



N 



(5(x 



X 



(7) 



Assuming that the statistics of the effective potential remains Gaussian, one finds within this first 
order perturbation theory: 



i?R(y) = Riy) + 



dl 

8^ 



didjRdidjR - didjRdidjR{0) 



(8) 



in D = 4, where b = e and dl is infinitesimal. 

Eq. (^) is the final result of the mode elimination. The search for fixed points requires the additional 
step of a rescaling transformation, which restores the original value of the cutoff A. Performing this 
rescaling via x ^ 6x and (j) ^ b'^cj) results in the full RG equation for the correlator 



diR = {e-4C)R + C(|>''^^R + 



1 

8^ 



didjRdidjR - didjRdidjR{0) 



(9) 



Iteration of this equation from the 'initial' condition R{y) = Ro{y) converges towards the fixed point 
R*{y), describing the long wavelength properties, which has the singular small y expansion [|lO| 



i?*(y)-i?*(0) =ey[ai-a3/2^] + ..., 



(10) 



where e = A — D \s the small parameter justifying the use of perturbation theory. Another way of 
stating this result is in term of the effective force f acting on the manifold, defined as minus the 
derivative of the effective potential with respect to (j). The force correlation function then behaves as 



[f*{<^)-f*m' 



12ea 



■3/2 1 



(11) 



Together with the assumption of Gaussian statistics, this suggests that the effective force acting on 
the manifold behaves, for = 1, as a random walk in space. This picture was advocated in [|lO| ], 

3 

and was actually used to argue that the next correction in e to would be of order ea . 



3 The replica approach. 

The replica approach is, in some sense, more ambitious, since it provides an explicit probabilistic 
construction of the effective disordered potential seen by the manifold. On the other hand, the method 
can only be controlled in the A — > oo limit, where a Gaussian variational Hamiltonian becomes exact 
[^ ]. Let us however stress right away that a Gaussian Hamiltonian in replica space does not mean 
that the actual effective potential which we wish to characterize has Gaussian statistics. As we shall 
indeed show below, this is not at all the case. 

Let us sketch first how the correlation function R{y) can be calculated with replicas and compared 
with the FRG. (More details can be found in |l2|, |8|). The average free-energy F = —^InZ = 

— -^In / P(/)exp[— is computed as usual as the 'zero replica' limit InZ = lim„_^o ^ ~^ ■ The average 
of Z'^ can be seen as the partition function of the following n— replica Hamiltonian: 

a=l \ / a,b 



<^(x))^ 



2AA2 



(12) 
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where an effective attraction between replicas has emerged from the disorder average. The idea is to 



treat this interaction using the trial Hamiltonian |26] 



^ ^EE r{-^)G-Jik)^{k), (13) 

a,b k 

where (/^"(k) = L~~ J d^:K (f)"'{x.)e~^^'^, and L is the 'linear' size of the manifold. 

The trial free-energy obtained with Hv depends on Gab and reads ^i,[G] =< Hn >v —-^TrlnG; 
the optimal matrix G is then determined by minimizing J^^ [G] , which leads to a set of self-consistent 
equations for Gab- The point now is that the structure of Gab in replica space can be non trivial 
in the limit n ^ 0, corresponding to 'replica symmetry breaking'. The physical meaning of this 
procedure has already been described in detail in [0, and we shall come back to it later. Before 

describing the solution to these self-consistent equations in the regime ~ 4, one should clarify first 
in what sense the replica calculation allows one to characterise the large scale pinning potential. Since 
the trial Hamiltonian is factorized over Fourier modes, one can isolate a particular, very slow mode 
ko — > 0. The effective force acting on ipQ = (^(kg) is f^{(pQ) = — ^^Vni^o), where Vni^Po) is the 
probability to observe (po for a given realisation of the random pinning potential It is thus clear 
that in order to compute, say, the correlation function of /, one should study the object: 



The last quantity is directly calculable, since the Gaussian Ansatz asserts that 



n -Pnm = j:^M-^0l^''>G'J{koWf\ (15) 



a=l 



where G is the optimal matrix determined via the self-consistent equations and vr denotes all the 
permutations of the replica indices. (All the saddle points only differing by permutation of the indices 
must be taken into account). The quantity in the right hand side of Eq. ([l^ ) corresponds to the 
choice (fQ = ifQ for ^ indices, and ipfj = ipQ for the other ^. The next trick to compute ( [T5| ) is to notice 
that in this case one can write (pQ = ^ [¥^0(1 + o'a) + Vo(l ~ ^a)], where Ca = il are fictitious Ising 
spins which pick up a particular permutation, provided J2a=i = 0. The technique for working out 
the sums over such spin configurations has been developped in the appendix (D) of ref.|^. Within a 
Parisi ansatz for the matrix G, the final result for the force correlation, written in the case of = 1 
to keep notations simple, is the following: 

2 A d'^ 1"^ 

fMhi^'o) = -,ckl-—j—— dh^{h,u = 0) (16) 



where ^{h,u) satisfies a non linear partial differential equation: 

where < u < 1 is the Parisi variable, indexing the pairs of replica indices a ^ b in the limit n ^ 0. 
The function q{u) is related to the matrix G(k, u) through: 

g(n) = -/3 ^^°"/°^' G-^(ko,n) (18) 



5 



and the boundary condition is 



$(/i,n = 1) = ln(^l + exp[-2/i-2g(l) + 2^ dn g(n)]j . (19) 

Hence, once G(k, u) is determined, the correlation function of the effective potential acting on mode 
ko is determined by solving (p!7|), which depends on (pQ — ip'^ through q{u). 

The solution of the self-consistent equations for G(k,u) was discussed in |11, Let us specialize 
to the case = oo, and introduce two important physical quantities, namely: 

• The Larkin-Ovchinnikov length (^lq separating a 'weakly distorted' regime for |x| < ^loj where 
all the displacements induced by the random potential are small compared to the correlation length 



of the potential A, from a strongly distorted regime. Simple dimensional arguments lead to pq, 12 



1 

„2 A 4 \ 4^ 



where is a rescaled potential strength, defined as = {2-k)'^W/{4: — D). The reason for introducing 
this rescaling in comes from the non trivial phase diagram around dimension D = 4 |29|. Indeed, 



it is easily seen from the study of the linearised, random force problem, that a 'weak disorder' regime 
with non trivial wandering exponent only exists when W is small enough. If one keeps the original 
W fixed and lets the dimension D go to 4, one enters a different phase (which actually survives for 



D > 4) |9|. 



• A 'Reynolds' number Re (this terminology comes from the analogy with Burgers' equation [Q), 
defined as the ratio of the elastic energy stored in a volume to the temperature ^. We shall define 
Re as 

Re = pC{D)wi=^ {cA^)^ , (21) 

2 

where C{D) is a dimension dependent number. Note that for D = A — e with e small, C = y- 

We shall only consider the case of low temperature and weak disorder, so that Re ^ 1 and ^lq ^ o, 
where a = ^ is the small scale lattice constant which regularizes the integrals over k. Under these 
conditions, we obtain the following result for Z) > 2: 

4-D 

G-\k,u) = -^1 ( — ) for u<Uc (22) 

= — (Ti for Uc < u < 1 (23) 

with (Ti = and Uc = 

The non trivial dependence of G^^(k, u) on u corresponds to continuous 'replica symmetry break- 
ing'. Let us now analyze the partial differential equation ( |l7|) in the limit /? — > 00. To this aim, we 
introduce the notation 7 = and the following rescaled variables: 

u = UcV il) = Uc^ g = I3aiul (^^^-^y) ('^o - 'f'of z = Uch/ (24) 
Eq. ( [l7|) then transforms into 

^ = -v''[i^" + vij'\ (25) 
(the ' means ^), with boundary condition (in the limit /3 — > 00): 

^{z, v = l) = -2 {{D - 2)g + z^) Q{-{D- 2)g - z^) , (26) 
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where is the step function. The correlation of free energies ( |16D thus involves, after change of 
variables, the integral 



fa 

dh^(h,u = 0) = 

uj. 



dz ziIj' {z, V = 0). 



(27) 



Under this form, the problem of evaluating I for small |(^o ~ '/^l can simply be treated by solving Eq. 
(25) for ip'(z,v), perturbatively in g - see Appendix A. The result reads: 



1 



7 + 1 



(28) 



with £(7) a complicated function of 7. In the limit D = 4 — e, 7 ~ — 1 + and £(7) — 2-^/7r. 
Transforming back to the original variables, we find, in the limit fco — > : 



RRSBiy) - Rrsb{0) = -eW 



y 



2A^ 



1 



^/y 



AC 



D/2 
LO ' 



(29) 



with y = {(po — (^o)^- Quite remarkably, Eq. (29) has the same form as the ERG result, Eq. ([lO|), 

— 3 — 

provided W is chosen in such a way that ^lq remains fixed as -D — > 4. This behaviour was first 

obtained within a replica theory in Q in the case D = 1 (corresponding to Burgers' turbulence), where 

the solution has a simpler, 'one-step' structure (valid for D < 2): G^^{k,u) = —aiQ{u — Uc)- 



4 Physical interpretation 

4.1 Shocks and relationship with the Burgers' equation 

As mentioned above, the Gaussian variational ansatz does not mean that the statistics of is 
Gaussian. Let us first discuss the replica construction of the effective potential in the simpler case 
D = 1 where a one step solution holds |11, ^. In this case, one has: 



y^*(^) = -iin 



(30) 



where a label the 'states', centered around ipa and of free-energy Fa, both depending on the 'sample' 
Q. The major prediction of the replica theory is that the Fa are exponentially distributed for 'deep' 
states [|, i.e: 

p{Fa) ccF^-^-oo eyi-p{-f3uc\Fa\). (31) 

The full distribution of the effective force (corresponding to the velocity in the Burgers problem) 
was analyzed in detail in [^]. Using the turbulence language, it was found that the velocity field 
organizes in a 'froth-like' structure of N— 1 dimensional shocks of vanishing width in the limit Re — > 00. 
Correspondingly, the potential has for = 1 the shape drawn in Eig. 1: it is made of parabolas 
matching at angular points - the shocks. The singular behaviour of the force-force correlation function. 



Eq. (11), is due to the fact that with a probability proportional to the 'distance' \(pQ — (^q|, there is a 



shock which gives a finite contribution to /((^o) ~ fi'^o)- This means in particular that all the moments 
\f{0o) — f{0o)\^ grow as \ipo — (^o| for p > I, instead of — as for Gaussian statistics. It is not 
clear how this strong departure from Gaussian statistics can be incorporated in an ERG treatment 
(see Section 4.C below). 

^This can actually be understood within the general context of extreme event statistics. [M 
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The relation with Burgers' equation is not coincidental and actually quite interesting. Keeping 
= 1 for simplicity, consider a toy model for the FRG mode elimination in which the renormalized 
effective potential is defined as 



/?^i?('/'<) =-ln 



(32) 



This means that Vr(v<) is precisely the Cole-Hopf solution of the Burgers equation |^l|: 

(33) 



di ~ 2/3cA2 9^2 ^ 
with 

V{ip, t = 0) = Vo{^) Vniip) = V{ip, t = 1). (34) 

As is well known a random set of initial conditions (here the bare pinning potential acting 

on (f) develops shocks which separates as time grows, between which the 'potential' V{ip) has a 
parabolic shape. Elimination of fast modes in a disordered system thus naturally generates a 'scalloped' 
potential, with singular points (which are smoothed out at finite temperature or finite Re) separating 
potential wells - the famous 'states' appearing in the replica theory. Quite remarkably, this structure 



was anticipated in [33, EQ] using different arguments. 



4.2 Full RSB and multiscale effective potential 

In the case of continuous RSB, the effective potential is recursively constructed via a set of 'Matrioshka 

doll' Gaussians. It is schematically drawn in Fig 2 for the the transverse fluctuations (j){£) — (p{0). 

For each length scale i, one can define a characteristic value of the parameter u{£) which plays the 

role of Uc in Eq. ( |3l| ) and sets the scale of the energy fluctuations. u{i) is such that the diagonal 

part of G~^{kQ = ^), namely cfcg, is equal to the off diagonal part G~^{ko,u), which gives u{£) oc 
e 

I (^ = D — 2 is the 'energy' exponent in the case = oo, and is related to the small 

u power-law behaviour of G~^{ko,u)). The large scale structure of the effective potential is thus 
a succession of parabolas of depth oc £^ , but this envelope structure is decorated by hierarchically 
imbedded parabolas corresponding to all the smaller length scales, between I and £,lOi beyond which 
the shocks disappear, since one enters into the effectively replica symmetric random force regime. 
The important point however is that small scale shocks are much more numerous than large scale 
ones and completely dominate the small y behaviour of RRSBiu)'- see Fig 2. This explains why the 
above result ( |29| ) is independent of k^ and only reflects the structure of G~^{k,u) in the vicinity of 
Uc, corresponding to /c ~ On the other hand, quantities like [(t>{£) — 0(O)]2 are dominated by 

the region where u ~ u{kQ = ^), corresponding to large scale moves. More precisely, the main 
contribution to [4>i£)) — 4'{^)V comes from minima separated by a distance £^ which happen to be 
separated by an energy gap smaller than the temperature [0, 17|. This occurs with probability 
oc X {pu{ko)) (see Eq. {^). 

In other words, the effective potential calculated within the FRG procedure involves an extra step 
which we have not performed within the replica construction, which is a coarse graining of the (p 
variables. In the FRG calculation, one restricts to configurations which are such that cp is constant on 
scales £, and scales as £'' 0. The correct choice of then ensures that there are only a few shocks on 
the scale £. As we now discuss in a rather conjectural way, this is perhaps why the FRG can still be 
controlled, the departure from Gaussian statistics being in some sense 'weak'. 

^ Removing all the modes fc > fco in the replica calculation leads to a correlation function Rrsb (y) indeed dominated 
by the vicinity of it(fco) oc feg. 
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4.3 The FRG in the presence of shocks 

To understand the emergence of shocks in the FRG picture, and to assess their impact on the pertur- 
bative procedure, it is useful to study the above toy model for the renormalization group, defined by 



Eq. (|3^), which amounts to discarding the internal degrees of freedom. Following Ref. |1C], we write 



Eq. (H) at zero temperature (and after a rescaling) as 

1 



Vk{<P<) = min -|</.>|^ + Vi<p^ + 0>) . (35) 



The validity of the perturbative minimization scheme was discussed in detail in Ref. ||10|], assuming 
Gaussian statistics for the random potential V. Errors occur in the perturbative minimization scheme 
due to an incorrect choice among multiple minima in the effective Hamiltonian for (p^ . For a Gaussian 
potential, there is an extremely dense set of such minima, and such an error occurs essentially with 
probability one. The FRG appears to be saved, however, because the magnitude of the resulting error 
in the energy is small (i.e. higher order in e). 

A rather different picture emerges if one assumes a smooth potential with shocks (i.e. slope 
discontinuities in V) spaced by 0(1) distances. To understand the limitations of the perturbative 
minimization scheme in this case, consider the extremal condition of the toy model, 

0> = -^'(0< + 0>). (36) 

In a scalloped (piecewise quadratic) potential, a perturbative solution in (/>> converges to the minimal 
energy in the local well containing (/>> = 0. For \ V\ small , this is indeed the global minimimum, unless a 
shock occurs within a distance |</'shock| < 0(|y|), as can be seen by examining the effective Hamiltonian 
for (py in the neighborhood of a cusp. Provided that a shock is present, however, the incorrect minima 
is chosen with a probability of 0(1), leading to a large error in Vji. Thus for the scalloped potential, 
instead of persistant small errors, the perturbative minimization scheme is typically correct, but suffers 
from catastrophic rare events that generate large errors with small probability. 

An interesting simplification occurs if one considers a periodic random potential V. Such periodic 
potentials occur in models of pinned charge density waves |^, and random anisotropy XY magnets |^] . 
It is straightforward to show that repeated applications of the toy model iteration drive the potential 
towards a form with a single symmetric cusp per period. ^ For such a symmetric form, the perturbative 
minimization scheme always converges to the correct (deepest) minima of the effective potential, i.e. 
the local minimum is always the global minimum. Within the toy model, then, the perturbative 
minimization scheme appears to be asymptotically exact. Although errors may accrue in early stages 
of the renormalization, these decrease as the length scale grows and the final fixed point form is exact 
- provided the perturbation theory is carried out to all orders, of course! That the FRG and replica 



methods lead to essentially the same results in this case was underlined in [13|. 

The FRG consists, as does any renormalization group, of two parts: the mode elimination (accom- 
plished via the perturbative minimization scheme) and the rescaling transformation. The toy model 
allows a detailed study, in a somewhat schematic way, of the former. Within this framework, the 
non-analyticity of R emerges in a natural way via the generation of Burgers' shocks. The toy model, 
however, completely neglects the internal degrees of freedom of the manifold, whose rescaling is crucial 
for the power-counting in the full FRG. In particular, this rescaling not only leads to the existence of 
a fixed point for R{4'), but also formally renders the higher cumulants of V strongly irrelevant. 

There appears to be a degree of competition between the mode elimination, which favors shocks 
and the corresponding highly non-Gaussian distribution for V, and the coarse graining and rescaling 

^This is because a potential with multiple cusps generically has multiple wells with non-uniformly separated and 
unequal minima, so that repeated minimization allows the surviving modes to remain only in the deepest well. 



9 



transformation, which tends to keep the density of shocks to a low value (at least for small e). A 
complete description, which is unfortunately not available to us at present, should properly balance 
these effects against one another. The special considerations applicable for the periodic potential 
discussed above suggest that the FRG may indeed be well controlled in that case. More generally, the 
full accommodation of shocks into the FRG remains a challenging open problem. 



4.4 The 1 + 1 Directed Polymer 

An explicit model where this construction actually does not require the use of replicas or of the FRG 
is the N = 1, D = 1 (Directed Polymer) case. From independent arguments [04| one knows that 
the effective potential Vx{4') acting on the 'head' of an infinitely long polymer (x oo) is a 'random 



walk' in (j) space: [la;(</') — ^xifp')]"^ 14' ~ ^'l- (Notice the difference with Eq. ([iTl), which concerns 
the force, and not the potential). In particular, there are no shocks in Vx{(j))- Shocks appear when one 
coarse-grains the description on a scale 6. Let us define a coarse-grained potential on an infinitesimal 
scale r] as 

Vx\p) ^ -^In (37) 

where /C is an arbitrary local 'filter'. Iterating this procedure a large number of times | produces an 
effective potential which, again, satisfies a Burgers equation, but now with a long range correlated 
'initial condition' V^.(0). As is well known shocks also appear in this case, with an average 

spacing growing as 5^/^. The distribution of distance d between shocks furthermore diverges for small 
d as p^, indicating that there are shocks on all scales smaller than 5'^^^. All these results can 



alternatively be obtained within the replica framework |25, 30 1 



5 Discussion and Perspectives 

We have shown in this paper that the FRG and RSB techniques are not contradictory but comple- 
mentary. They both suggest quite an appealing physical picture: the phase-space of the system is, on 
large length scales, divided into 'cells' corresponding to favourable configurations where the potential is 
locally parabolic, and whose depth is exponentially distributed. These cells are themselves subdivided 
into smaller cells, corresponding to larger length scales, etc.. This hierarchical construction is similar 



to the one usually advocated for the phase space of spin-glasses [35|, based on Parisi's RSB solution 
of the SK model [|ll|, . The enormous advantage of random manifolds is that this construction can 
be directly performed in physical space. 

An important consequence of this construction is that it allows us to discuss the dynamical prop- 
erties for finite |^. In the case of a one-step RSB, one can directly calculate from Eq. (^), the 
distribution of the height of the barriers AE between two neighbouring wells, and finds that it decays 
exponentially as exp(—(3ucAE). It is interesting to notice that the barriers thus behave in the same 



way as energy depths [^], a point recently studied in detail for randomly pinned lines in |38]. A 
natural picture for the dynamics is thus to imagine that the manifolds jumps from well to well, each of 
which representing a long-lived conformation of the manifold. Such a picture is corroborated by recent 



numerical simulations in D = 1, N = I |39]. The lifetime of each 'trap' is activated r ~ tq exp(/3A£/), 
and is thus distributed as a power-law 7--i-"('=) for large r, where the exponent u{k) tx depends 
on the 'size' of the jump (i.e. the mode involved in the change of conformation), small u{k) corre- 



sponding to large wavelengths. Then, as emphasized in |1£] where precisely the same 'trap' picture 



^Infinite A'' properties have been thoroughly investigated in jS^ but might belong to a different dynamical class since 
activated effects disappear, at least for short range correlations 
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was advocated for spin-glasses, the dynamics becomes non stationary and aging effects appear at low 
temperatures and/or long-wavelengths such that u{k) < 1. For example, the response of the manifold 
to a spatially modulated external field is expected to behave, for t <^ ty^, as (^)^~"('^)^ where tw is 
the time elapsed since the quench from high temperature. Correspondingly, the a.c. response should 
behave, for wt^ ^ 1, as (wt^;)"^^^"^, again much in the same way as observed in spin-glasses [p^]. 
For finite N however, one may expect that the exponential distribution of deep states ceases to be 
valid outside the scaling region, i.e. for AE >>> ^J^^^ p9| , |30[| . This will lead to 'interrupted aging' 

for modes such that Int^ >>> u{k)~^. These equilibrated modes thereafter only contribute to the 
stationnary part of the response (or correlation). 

It is thus rather satisfactory that the 'traps' appear naturally in the context of pinned manifolds 
through the replica description, and that this picture actually complement the 'droplet' construction. 
It would of course be gratifying to understand precisely how these ideas could be extended to finite 
dimensional spin-glasses. 
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A Perturbation expansion for Eq. (27) 

We provide here some intermediate steps of the computation of free energy correlations with the replica 
method. We need to solve eq.(p5|) with the boundary condition (p^), and compute then the integral I 
defined in (|27|). The limit of interest is g small. We work with the derivative xi^^^) = '^'{z-.v) which 
satisfies the equation: 

^ = -v'[x" + vx^] (38) 

together with the boundary condition: 

X(z, v = l) = -2^e {-{D - 2)g - z^) . (39) 
The solution to this differential equation to order g can be written as: 

X{z,v = {))=xi{z)+X2{z) (40) 



where 



and 



(Z - 



Xi{z) = -2^g I -j^ exp - (41) 

^ \ 7+1 / 



/■I dz' {z-z'f \ d / fO dzr / {z - z'f 



7+1 \ 7+1 / \ V -y+l \ 7+1 

(42) 

Introducing the notation 

Moix) = f -^e-"^/^ (43) 
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we find, after multiplication of x by z and integration: 

(44) 

Expansion of the last integral for e = 4 — £) small, with 7 = — 1 + e/2, reveals that the coefficient of 
5^/^, which to leading order should be oc e~^^^ in fact vanishes, the next term being of order e°. 

Figure Captions 

Fig. 1. Schematic view of the effective energy landscape as a succession of parabolic wells matching 
at singular point. This picture actually corresponds to a 'one-step' replica symmetry breaking scheme. 
Fig 2. Multiscale energy landscape corresponding to a full replica symmetry breaking scheme. In 

this case, the construction is that of parabolas within parabolas, in a hierarchical manner. The depth 

e 

of the wells (and thus also the height of the barriers) typically grows as \(f> — (/) \'^ . The figure actually 
corresponds to a two-step breaking scheme, with ui = 0.5 and uq = 0.05. The inset is a zoom on a 
particular region, showing the first level of Gaussians. 
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